Multi-electron transitions induced by neutron impact on helium 
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We explore excitation and ionization by neutron impact as a novel tool for the investigation of 
electron-electron correlations in helium. We present single and double ionization spectra calculated 
in accurate numerical ab-initio simulations for incoming neutrons with kinetic energies of up to 
150 keV. The resulting electron spectra are found to be fundamentally different from photoioniza- 
tion or charged particle impact due to the intrinsic many-body character of the interaction. In 
particular, doubly excited resonances that are strongly suppressed in electron or photon impact 
become prominent. The ratio of double to single ionization is found to differ significantly from those 
of photon and charged particle impact. 
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Spectroscopic studies of atoms, molecules, and solids 
rely on the well established excitation processes such 
as photoabsorption and charged-particle impact. The 
underlying dynamical processes are theoretically well- 
understood within the framework of linear response of 
the system to the external probe. The observables acces- 
sible by these probes are, however, limited by either ex- 
act selection rules or approximate "propensity" rules. For 
example, photoabsorption spectroscopy is strongly dom- 
inated by dipole-allowed transitions. In charged-particle 
impact, higher multipole transitions are allowed but are 
typically suppressed in "soft" collisions with small mo- 
mentum transfers. Moreover, the long-range Coulomb 
interactions between the probing particle and the ex- 
cited system may distort the excitation and ionization 
to be extracted by "post-collision" interactions which are 
typically beyond lowest-order perturbation (LOP) theory 
underlying linear response. 

Photon and charged particle interactions have in com- 
mon that the LOP interaction is strictly a one-body oper- 
ator. The point of departure of our present study is the 
observation that neutron impact gives rise to intrinsic 
many-body interactions in the electronic system [1]. The 
underlying idea is that neutron scattering at the atomic 
nucleus gives rise to a sudden "kick". In the frame of the 
atom, this results in a simultaneous momentum boost for 
all electrons, effectively causing a true many-body transi- 
tion which can efficiently lead to multiple excitation and 
ionization of the atom. In this Letter we theoretically in- 
vestigate the neutron-impact ionization of helium atoms. 
Helium is the prototypical case of a strongly correlated 
system [2] in both the ground state and doubly excited 
resonances which can be treated exactly by numerical ab 
initio calculations (cf. e.g. [3-6]). It thus serves as testing 
ground for the study of electron correlation and multi- 
electron effects. We show that neutron impact leads to 
a very broad energy distribution in the final states in- 
cluding double ionization and, furthermore, that it can 
efficiently produce doubly excited states that are disfa- 



vored by other probing agents. 

We assume that the only interaction in the neutron- 
helium collision is quasi-elastic scattering between the 
neutron and the nucleus, mediated by the strong nuclear 
force. The contributions of magnetic interactions of the 
neutron with the electronic and nuclear magnetic mo- 
ment are small enough to be safely neglected [7] . Neutron 
energies are kept sufficiently low in order to exclude any 
inelastic nuclear processes. The duration of the neutron- 
nucleus scattering event is much shorter than the typical 
time scale of the dynamics of electrons bound to the nu- 
cleus (~ attoseconds) currently probed using ultrashort 
light pulses [8-10]. Electronic transitions can therefore be 
described by an impulse or "sudden" approximation. Ac- 
cordingly, the transition amplitude for quasi-elastic scat- 
tering of the neutron accompanied by an electronic tran- 
sition i — > f is given by 



t e n UC (&Pnuc)-tl f (Ap e ) 



(1) 



where t„ uc is the transition amplitude for elastic nuclear 
scattering with momentum transfer Ap nuc — kf — ki and 
t\ j is the matrix element of the collective boost operator 



tl f (Ap e ) = (V f \ exp[iAp e ■ (n + r 2 )]|^, 



(2) 
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FIG. 1. Projected two-electron momentum distribution 
P(ki :Z , k2, z ) for (a) the ground state of helium, (b) the ground 
state wave function boosted by the one-body boost operator 
Bib with a momentum transfer Ap c = 1.0 a.u., (c) boosted 
by the collective boost operator B c (Eq. 4) with identical Ap e 
(see text). 
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with 



(3) 



M a +2 

and M„ the mass of the a particle in atomic units. Taylor 
expansion of the collective boost operator 



B c (Ap e ) = exp[iAp e ■ (f x + r 2 )\ 



1 + iAp e • (rl [Ap e • (ri + f 2 )] 2 



(4) 
i 

(5) 



shows that while, to first order in Ap e , the electronic 
transition matrix element is equivalent to that of the 
one-body operator from photoabsorption or the Bethe- 
Born limit of soft charged-particle collisions, all higher- 
order terms represent a true many-body transition struc- 
turally different from photon or charged-particle interac- 
tions. Application of the collective boost to the exact he- 
lium ground state (Fig. 1) leads to a correlated displace- 
ment of the projected two-electron momentum distribu- 
tion unlike the one-body boost operator, _BiB(Ap e ) — 
y~lj—i exp(iAp e • fi) governing, for example, Compton 
scattering or charged particle impact on an A^-electron 
atom. This property plays a key role in accessing states 
blocked by parity or propensity rules. 

Differential cross sections for electronic inelastic pro- 
cesses accompanied by quasi-elastic neutron-alpha parti- 
cle scattering are given by 

dcr,^f , . . kf da, 



dQ 



l (Ap„ 



ki dfl 



(Ap nuc )\t^JAp 



(6) 



with kf = y/kf — 2fiQi. Qi = Ej — Ef the internal exci- 
tation energy, and fj, the reduced mass of the n-He system. 



For the nuclear elastic scattering cross section 



do 



use the tabulated data from [11]. For the electronic de- 
grees of freedom in helium we perform full ab-initio calcu- 
lations by solving the six-dimensional time-independent 
Schrodinger equation (five-dimensional after exploiting 
cylindrical symmetry) including all interparticle interac- 
tions. In our computational approach we employ a close- 
coupling scheme, in which the angular variables are ex- 
panded in coupled spherical harmonics (with total angu- 
lar momentum up to L max = 7, and individual electron 
angular momenta up to Z max = 9). For the discretiza- 
tion of the radial components we use a finite element 
discrete variable representation (FEDVR) [12, 13]. The 
momentum boost operator Eq. (5) is implemented us- 
ing a short iterative Lanczos algorithm (SIL) [14]. For 
the extraction of transition amplitudes, the direct pro- 
jection onto final states would be most desirable but un- 
feasible as exact three-body Coulomb continuum states 
are not known. We therefore make use of an alternative 
approach [5, 15], in which the Fourier transform of the 
boosted wave packet is effectively calculated by solving 
the inhomogeneous linear system 
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FIG. 2. (Color online) Ratio of double to single ionization 
for neutron-impact ionization Rn(Ae) as a function of en- 
ergy transfer Ae compared to the corresponding ratio Rc for 
Compton scattering [16] multiplied by a factor 6 for visibility. 
Inset: The same ratios normalized to their respective maxima. 
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(E - H)\V SC (E)) =B c (Ap e )\y z ), 
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FIG. 3. (Color online) Single (SI) and double ionization (DI) 
probabilities of helium as a function of the momentum boost 
Ap e for the electrons. The corresponding recoil energy of the 
kicked helium nucleus is given on the upper abscissa. 



where ^> SC (E) is the scattered wave function in the (time- 
independent) energy domain. Outgoing boundary condi- 
tions are enforced by an exterior complex scaling (ECS) 
transformation for each of the radial coordinates. For the 
calculations presented in this Letter we chose an exterior 
scaling radius of 120 a.u. and an overall box size of up to 
180 a.u. . The ejected single and double ionization ampli- 
tudes can then be extracted from the scattering ampli- 
tude by means of a surface integral within the non-scaled 
part of the grid [15]. 

The most frequently studied quantity in double ioniza- 
tion of helium, a paradigm for studying the role of elec- 
tron correlation, is the ratio of double to single ionization 
R. This ratio has been probed for both charged particle 
impact and photon impact over a wide range of energies, 
both experimentally and theoretically [17]. For photon 
impact, photoabsorption as well as Compton scattering 
have been studied [16, 18-23]. Compton scattering in- 
volving a neutral projectile and the one-body boost op- 
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FIG. 4. (Color online) Absolute integrated double (odi) and 
single ionization (<tsi) cross section by neutron impact as a 
function of the neutron kinetic energy. Inset: Energy de- 
pendence of _Rn = o"di/o"si with magnification of threshold 
region. The nonrelativistic high-energy limits for photoab- 
sorption 7?pa, Compton scattering _Rc, and charged particle 
impact Rc are shown for comparison. 



erator as transition operator is expected to bear closest 
resemblance to the present case of neutrons. Significant 
differences are, however, expected, as for neutrons the 
collective boost rather than the one-body boost controls 
the transition and, moreover, different regions in the en- 
ergy transfer (Ae) - momentum transfer (Ap e ) plane are 
sampled. The ratio for neutrons, R^(Ae), differential in 
energy transfer Ae to the electronic system, qualitatively 
resembles the calculated Rc(Ae) for Compton scattering 
near threshold (inset Fig. 2). Its absolute magnitude is, 
however, strongly enhanced by factors up to 25 depend- 
ing on the kinetic energy of the incident neutron (Fig. 2). 
The most dramatic difference (Fig. 3) occurs for large 
momentum transfers due to the non-linear dependence 
of the boost operator on Ap e . In the limit Ap e — > oo, 
or more precisely when the momentum transfer is large 
compared to the width of the momentum distribution of 
the initial state, Ap e ^> (pf) 5 , the ratio diverges, as the 
strongly displaced momentum distribution (Fig. 1) will 
effectively cease to overlap with bound states and double 
ionization dominates. 

Most easily accessible in experimental investigations is 
the ratio i?N of total double to single ionization cross sec- 
tion (Inset Fig. 4) resulting from integration of Eq. (6) 
over all accessible final states in the Ae-Ap e plane as 
a function of the kinetic energy of the incident neu- 
tron. With increasing neutron energy the ratio Rn in- 
creases polynomially (cx aiE^ + a 2 -E^ + • • • ) with the 
neutron energy and eventually surpasses the well known 
(non-relativistic) high-energy limits for photoabsorption 
(1.66%) [18-21], Compton scattering (0.8%) [16, 22, 23] 
and charged-particle impact (0.26%) [17, 24]. The rea- 
son is that the He nucleus suddenly "disappears" from the 
electronic charge cloud resulting in a high probability for 
double ionization. 
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FIG. 5. (Color online) (a) Electron spectrum as a function of 
final energy for single and double ionization for a kick strength 
of Ap e — 1 a.u.. The energy of the ejected electron in the case 
of SI is given by E + \ with the remaining ionized helium 
being excited to state n. In the case of double ionization 
the probability for ejecting electrons with the sum of the in- 
dividual electron energies equivalent to E is plotted. The 
full two-electron energy distribution is plotted in the inset 
of (a). A close-up of Fano resonances in the n — 1 channel 
for different final symmetries is shown in (c). The first few 
resonances are labeled by approximate independent-particle 
configurations. For comparison, the one-photon spectrum (b, 
calculated) and an ejected electron spectrum (d) from e~-He 
scattering experiments performed by deHarak et al. [25] are 
shown. 



A more sensitive probe of the momentum shift of the 
two-electron momentum distribution by the collective 
boost is the energy spectrum of the ejected electrons 
in single and double ionization (Fig. 5). The impul- 
sive momentum transfer leads to a broad-band excita- 
tion (the upper cut-off due to the finite nuclear collision 
time ~ 1/tcoii lies well beyond the spectral range shown 
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in Fig. 5) resulting in a large number of doubly excited 
resonances embedded in the single ionization continuum. 
A zoom into the electron energy spectrum just below 
the n = 2 threshold [Fig. 5(c)] shows the multitude of 
Beutler-Fano resonances of different symmetries. Note 
that doubly excited states are not properly identified by 
the usual independent-particle labels, but require collec- 
tive quantum numbers (cf. [2] and references therein). 
However, for brevity, we use the traditional but impre- 
cise labels (nln'l') to describe the first few doubly ex- 
cited states. The background from direct single ioniza- 
tion into the continuum is only strong in the channel with 
1 P° final symmetry, while it is suppressed in the other 
channels. This is a clear signature of the different domi- 
nant terms in the transition operator for different symme- 
tries: in 1 P°, the first-order (one-body) part dominates, 
which couples efficiently to the single continuum, but 
only weakly to doubly excited states. In 1 S e and 1 D e , the 
dominant part of the boost operator is the second-order 
two-body term. The latter couples the initial ground 
state more efficiently to the quasi-bound doubly excited 
states than to the single ionization continuum. This is 
best seen in the (2p) 2 (both in the 1 S' e as well as in 
the 1 D e channel) and (2p3p) doubly excited states which 
feature the largest cross section. By contrast, these tran- 
sitions are strongly forbidden in photoabsorption driven 
by the dipole operator (first term in Eq. 5). Exciting 
those resonances by photons would require a two-photon 
absorption process triggered by an intense beam with 
well-tuned frequencies, in reach with free-electron lasers 
[26, 27]. Even within the dipole-allowed 1 P° spectrum 
neutron-impact ionization leads to a marked modification 
of the Beutler-Fano resonance profiles [28, 29] compared 
to photoabsorption [Fig. 5(b), (c)]. The latter is a sig- 
nature of interference between the first and third-order 
terms in Eq. 5. 

It is instructive to compare the neutron-impact in- 
duced spectrum with the corresponding spectrum for 
electron impact [Fig. 5(d), taken from [25]]. While for 
charged-particle collisions higher multipole transitions 
become allowed, the propensity for excitation of res- 
onances of different symmetry arc markedly different. 
Considering, for example, the first two doubly excited 
resonances in the 1 S e channel, the (2p) 2 state is much 
stronger excited for neutron impact ionization than for 
electron scattering. This is in contrast to the (2s) 2 dou- 
bly excited state, which is present in both excitation pro- 
cesses. This difference can be explained by specific elec- 
tron correlation effects present in these doubly excited 
states. It has been shown [30] that a major difference 
between the two states lies in the expectation value of 
the angle #12 between the two electrons. For the (2s) 2 
state the electrons are more likely situated opposite to 
each other whereas in the (2p) 2 case they have a ten- 
dency to be located on the same side of the nucleus. For 
the quasi-instantaneous neutron kick it is suggestive that 



both electrons will be pushed to the same side of the 
nucleus and will thus have significant overlap with this 
class of resonances. This behavior is less likely for exci- 
tation by an incoming electron which interacts with the 
bound electrons via the long-ranged Coulomb force and 
gives rise to transition matrix elements containing the 
one-body boost operator. The strong excitation can thus 
be directly attributed to the effective many-body nature 
of the neutron kick. In contrast to neutron impact, the 
collision with an incoming electron can also access 3 P° 
states due to spin exchange processes, which can be seen 
in Fig. 5(d) for the 3 P°(2s2p) state. 

Analogous processes induced by neutron impact are of 
interest also in larger systems, e.g. in molecules, where 
they result in the excitation of auto-detaching states and 
opening of dissociative channels. In solids, they represent 
the key processes underlying electronically induced radi- 
ation damage triggered by energetic neutrons subsequent 
to knocking nuclei from their lattice positions. 

In conclusion, we have shown that neutron-impact 
ionization could serve as a novel tool to probe corre- 
lated electronic dynamics in many-body electron sys- 
tems, specifically in helium. Key is the true many-body 
nature of the correlated boost operator which allows tran- 
sitions that are either strictly forbidden or strongly sup- 
pressed in either photoabsorption or charged-particle ex- 
citation. Doubly excited resonances become prominent 
that are otherwise only barely visible. The ratio of double 
to single ionization by neutrons, R^. is another bench- 
mark for the underlying differences of the ionization pro- 
cess. The predicted ratios significantly differ from those 
for photoabsorption, Compton scattering, and charged- 
particle collisions. With the availability of high-intensity 
neutron sources, the observation of these processes under 
well-characterized single-collision conditions may come 
into reach. 
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